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ABSTRACT 

The equation of state of symmetric nuclear matter with the inclusion 
of non-strange dibaryons is studied. We pay special attention to the exis- 
tence of a dibaryon condensate at zero temperature. These calculations have 
been performed in an extended quark-meson coupling model with density- 
dependent parameters, which takes into account the finite size of nucleons 
and dibaryons. A first-order phase-transition to pure dibaryon matter has 
been found. The corresponding critical density is strongly dependent on the 
value of the dibaryon mass. The density behavior of the nucleon and dibaryon 
effective masses and confining volumes have also been discussed. 

PACS: 21.65.+f; 12.39.Ba; 24.85.+p; 14.20Pt 

Keywords: Nuclear matter; Bag model; Dibaryons; Equation of state; Bose condensate. 



*Electronic address: aguirre@venus. fisica.unlp.edu. ar 
'I' On leave from University of La Plata. 
Electronic addresses: schvell@cebaf.gov and martin@venus.fisica.unlp.edu.ar 



Nuclear matter properties have been modeled by several kind of local relativis- 
tic effective lagrangians which basically use point-like representations of nucleons and 
mesons as the relevant degrees of freedom. The differences are essentially coming from 
the nucleon-scalar meson interaction. It should be desirable to describe nuclear matter 
from a fundamental theory like QCD, but it is wellkown that in the low energy limit 
QCD becomes non-perturvative. This situation has motivated the development of several 
effective models of quark interaction, among them there is the MIT bag model [0. This 
model predicts the existence of some multibaryons and strange exotics 0]. The early 
work of Jaffe has concentrated the attention of theoretical studies [0, ^ and the fur- 
ther development of experimental research looking for signals of strange and non-strange 
dibaryons. In reference it has been shown that in the scheme of the MIT bag model 
the gluon-exchange force should be responsible of the existence of a stable six-quarks 
bound state. This particle, the so-called if-particle, is a flavor singlet J'^ = dihyperon 
with a mass of 2150Me\^ and strangeness —2. The problem has also been treated in the 
non-perturbative framework of the Skyrme model where dibaryons have been considered 
as axially symmetric skyrmions |^ . Non-strange as well as strange dibaryons and another 
multibaryons have been studied as multiskyrmions 0. The experimental activity con- 
cerned to this search has been increased in the last years, recent experiments have been 
developed at TRIUMF and CELCIUS 0. Non-strange dibaryons which have a small 
width have been described as promising candidates for experimental searches 0. 

On the other hand, models based on the quantum field theory of hadrons including 
non-strange dibaryons as effective degrees of freedom have been extensively studied P], 
obtaining very interesting effects on the binding energy per particle as well as on the 
equation of state (EOS) of the system. In these hadronic effective lagrangians nucleons 
and dibaryons are treated as point-like particles represented by two independent effective 
fields, both of them interact by the exchange of scalar and vector neutral mesons using 
two different sets of coupling constants. 

The purpose of this work is to include finite size effects on the EOS, taking into 
account the quark structure of the particles. As we will show later, several features of our 
results are in agreement with those obtained in reference 0. Therefore, this fact seems 
to support that the dibaryon condensate is essentially a model-independent issue. We 
have selected the so-called quark-meson coupling model (QMC) [|I^, which in some way 
satisfies the above mentioned requirements. In this scheme we can perform a simultaneous 
description of nucleons and dibaryons using only current quarks and mesons as effective 



degrees of freedom. The model was early developed by Guichon |10] and it has been 
extensively applied to calculate nuclear matter |]ll| as well as finite nuclei [|l^ properties 
with successful results. It has also been used to evaluate nucleon structure functions 
1131 . The naturalness of the QMC model has been studied using the dimensional analysis 



Tj]. Recently the density dependence of the parameters of the MIT bag have been 
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phenomenologically modeled |T^, alternatively a relationship with observables evaluated 
in the quantum field theory of hadrons have been stablished . 

In this work we have developed an extension of the QMC model described in refer- 
ence [|l6l, by including dibaryons represented as spherical MIT bags confining six quarks. 
The extended QMC lagrangian density with quark fields qa{x) coupled to scalar a{x) and 
vector uj^{x) neutral mesons, is written as follows 

CqMc{x) = Cn{x) + Cd{x) + C%^sonsix) , (1) 

C^ix) = {C%{x) - B,) QvN - i ELi qa{x)q^{x)AsN , (2) 
Coix) = (CUx) - B2) QvD - \ ELi ga(a;)g„(x) A^d . (3) 

Here and Oyz? are the radial non-overlapping step functions which schematically 

confine the quarks inside spherical bags for nucleons and dibaryons, respectively. B\ 
and B2 are the so-called MIT bag constants associated with these particles. Within the 
standard QMC treatment S is a constant, however it can be considered as function of 
the baryonic density, ps- The terms proportional to the surface delta functions A^^v and 
l^SD ensure a zero fiux of quark current through the bag surface. In this lagrangian we 
have defined the following terms 

C^%{x) = ELi Qaix) {iYd^ -771^+ gaf^ix) - g^Y^f^ix)) Qaix) , (4) 

^oi^) = ELi Qaix) {i^d^ -ma + ga<y{x) - g^'^^ujf,{x)) Qaix) , (5) 
^Mesonsi^) = '^[d^cr{x)d,a{x) - mla\x)] - \F'^^{x)F^,{x) + \mluj^{x)uj,{x) . (6) 

Here g^, and gi^ are the quark-meson coupling constants associated with ct{x) and co'^(x), 
respectively. In order to get the minimal set of free parameters we have assumed that 
quarks do not distinguish baryon or dibaryon bags, i.e. the coupling constants are the 
same in both of the equations for baryons and dibaryons. In that follows we deal with u 
and d massless quarks, furthermore we introduce the index v = 1, 2 to label quantities 
related to nucleons and dibaryons, respectively. 

The normalized quark wave function for the fundamental state in a spherical bag 
of radius Ry is given by 

g.(f, t) = ACe-*/-^ X ( '±'1''^]^ (7) 

where r is the distance from the center of the bag, Xv is the quark spinor and the nor- 
malization constant is 

K = , . (8) 

^2Rlji{y,)[n,{n, - 1) + R,m*j2] 

The parameter associated with the quark mass is m* = — g^a, and the energy eigen- 
value is written as 

= ^lu + g^jVOjRu , (9) 
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where VL^ = ^y1 + (i?,^m*)2, and a,uj are the mean values of meson fields calculated in 
the Mean Field Approximation. The yy variable is fixed by the boundary condition at the 



bag surface jo(l/i/) = l3uji{yu) as in reference ^ and (3^ = \J {^u - Ryml)/{VLy + Rural). 
The mass associated with the bag is given by 

= + 4 5^^3. (10) 

The By are the bag constants already introduced in Eqs.(2) and (3), while zq^ takes into 
account the zero point energy of the bag. 

The usual procedure in QMC is to fix the nucleon bag parameters at zero baryon 
density to reproduce the experimental nucleon mass Mi = = 939MeV, simulta- 
neously it is required that the equilibrium condition dMi{a)/dR = must be fulfilled. 
An analogous method should be applied to dibaryons, however the experimental value of 
in- vacuum dibaryon mass has not been definitely confirmed yet, and at the present, 
only theoretical estimates are availables. Therefore, we leave as a parameter of the 
model. Since the B^, are related to vacuum properties we assume that Bi = B2 = B, 
at all densities. Although the paramaters Zq^ could also have a density dependence, in a 
previous work [16] it was found that zqi remains approximately constant at the baryon 
densities below four times the nuclear matter saturation density. Consecuently we have 
taken zq^, as constants fixed at zero density for each kind of bag. Under these assumptions 
one can immediately get a relation for masses and radii : R2 = {M2/Mif/^Ri. 

In order to obtain the density dependence of B we have stablished an explicit 
relationship between nuclear matter observables evaluated in the extended QMC model 



and in pure hadronic models, as it has been described in detail in [IC]. In the present work 



we have selected the Zimanyi-Moszkowski model [|18| to describe the hadronic sector (as 



in reference [|TB|). Considering that the effective nucleon mass predicted by the extended 
QMC model and the hadronic model must be the same, together with the fact that 
the outward quark-momentum on the bag surface must be compensated by the hadronic 
momentum going inside, one gets the following equations 

Ml = M* , (11) 

Pbag = Phad , (12) 

which are valids at each value of baryon density. Here M* and Phad are the effec- 
tive nucleon mass and pressure of nuclear matter in the hadronic model, while Pbag = 
— {dMi/dRi)/4:7TRl is the pressure in the nucleon bag. Using these equations and mo- = 
550MeV, = 783MeV, = 4.576, = 2.222 and P? = O.Sfm we have obtained the 
density behaviour of B = B{pb) shown in Fig. 

Once we have dinamically derived the parameter i? as a function of the baryon 
density, we can perform the calculations for nuclear matter with dibaryons. Firstly, we 
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evaluate the energy per baryon of a system composed by symmetric nuclear matter with 
density p^-, and dibaryons with density po- Since each dibaryon carries baryon number 
two, the baryon density must be ps = Pn + '^Pd- Using the quark wave functions of Eq.(|^ 
one can construct the antisymmetrized (symmetrized) nucleon (dibaryon) physical states 
and evaluate the expectation value of the energy density Ti derived from the lagrangian 
of Eq.(l). For uniform matter the energy per baryon is given by 



e = 'H/pB = — 

PB 



(13) 



where r] = kp/Mi, kp is the Fermi momentum for the nucleons and it is related to the 
nucleon density by pn = 2kp/3n'^, while F{ri) = rjy/r]'^ + l{2rf + 1) — Logij] + ^rj^ + 1). 



Chemical potentials for nucleons and dibaryons are given by p^ = ykp + Mf + 3gi^Lj 
and Pd = M2 + 6g^uj, respectively. Although the dibaryons are bosons, the corresponding 
particle number is conserved because they carry baryon charge and since in our lagrangian 
we have not included any decay mechanism. The pressure of hadronic matter at zero 
temperature can be written as 

Phad = PnPn + PdPd - epB ■ (14) 

The mean values of the meson fields have been evaluated by minimizing the energy, i.e. 
de/da = and de/duj = 0, obtaining 



cr 




dF 

AMiF{r]) - pF—{r]) 
arj 



dMi , dM2\ 

-^+PD-^\ , (15) 
da da 



together with lj = Sg^pB/rri^. 

In order to obtain numerical results we have arbitrarily selected the value of the 
in-vacuum dibaryon mass = 1970Mel^. We have used the meson masses, free nu- 
cleon bag radius and the previously described procedure to calculate several in-vacuum 
quantities such as the parameters 2:01 = 3.27, 2:02 = 6.27, as well as the dibaryon bag 
radius = 0.954 fm~^, the usual boundary condition eigenvalue yj^^ = 2.042 and the 
bag parameter = 0.5546 fm~'^. When the constants Zqi? ^02 have been fixed and the 
density dependent parameter B has been obtained, we must find the effective bag radii 
Ri{Pb) and R2{pb)- To do that we have impossed the equilibrium condition for bags 
immersed in hadronic medium in similar way as Eq. (p!2D , but now using Eq.(p!4D in the 
right hand side. After a little algebra it gives 

^^RtPkad = -SM^R. + 4(3fi. - .0.) - 3«.^^^^P^\^ , (16) 

where we have defined a^, = m'^Ri,{pB). It must be noticed that this condition becomes 
the standard equilibrium requirement of QMC model for vanishing pb- 
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Solving the boundary condition at the bag surface simultaneously with Eqs. ([15|) and 
(p!6|), for a given baryon density and several values of the dibaryon abundance Q = Pd/pb 
we have evaluated the bag radii, masses, energy density and pressure as functions of 
and Q. The coupling constants have been fixed in order to reproduce at Q = the 
saturation density po = 0.15 fm~^ and the binding energy per particle eo = 16MeV, 
obtaining = 4.6876 and = 2.2967. In Fig. 1 the results for the nucleon bag radius 
Ri are shown. From this we can see that at low densities the radius is an incresing function 
of pb and it becomes a smoothly decreasing function at densities around the saturation 
nuclear density. Furthermore, at densities below 1.4po, -Ri is higher than in- vacuum radius, 
notwithstanding for Q = the relative increment of -Ri remains under the value 2%, as 
it is expected from y-scaling arguments in the analysis of quasi-elastic electron scattering 
||19|| . Therefore, the use of this radius and the nucleon effective mass Mi is equivalent to 
have a bigger radius and the free nucleon mass, as it has been proposed in the reference 
of Sick [0. Density variation of the radius is less pronounced as Q increases. A similar 
behaviour for the dibaryon bag radius has been found. 

To investigate the existence of a phase-transition from symmetric nuclear matter 
to a pure dibaryonic state we have used the Gibbs criterion. If there is a dynamical 
mechanism combining two nucleons to give one dibaryon, at the phase-transition point 
the conditions 2pj^ = p^ and Pnuc = Pdib must be fulfilled. Although we have not 
included such a reaction channel in our model, we are concerned with the stable initial 
and final phases rather than with any particular mechanism of dibaryon formation. Here 
Pnuc and Pdib denote the pressure in pure nuclear and pure dibaryon matter, respectively. 
In Fig. the quantities 2pN for Q = and pn for Q = 0.5 are represented as functions of 
the corresponding pressures. The intersection point of these curves indicates that phase- 
transition occurs at a critical pressure nearly Pq = 8.8 10~^/m~^. At low pressure the local 
minimum of the Gibbs potential per particle corresponds to pure nuclear matter, while 
at pressures beyond Pq the stable state corresponds to pure dibaryon matter. In order to 
construct the EOS we have drawn Fig. ^ with the pressure of Eg . ([l^ , corresponding to 
Q = and Q = 0.5 as functions of baryon density. The value Pq is reached at baryon 
densities of p/ = 2.55po in nuclear matter and pu = 3.18po in the dibaryon condensate. 
Therefore it is a first-order phase-transition with a discontinuity jumping 0.63po in the 
baryon density. The horizontal segment P = Pq between p/ and pn represents two-phase 
coexistence. The effect of the phase-transition is to reduce the compressibility at high 
densities. For example, the quotient of the thermodinamical compressibility evaluated in 
the dibaryon condensate to its value in the nuclear matter gives at ps = 3.5po the value 
0.706. This fact could be favorable for the collapse of very massive stars [^. In Table I 
a comparison of the quotient of thermodinamical compressibility at Q = 0. and Q = 0.5 
for several values of Md at the transition point is shown. 

In addition, we have studied how much is modified the phase-transition point as 
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the in-vacuum dibaryon mass Md is changed. It has been found that the density pi is 
an increasing function of the mass M^, as it is shown in Table I. We must take into 
account that dibaryons have not been observed at densities around the normal saturation 
density po and at sufficiently high densities our model can not describe the quark-gluon 
plasma phase-transition. Therefore we limit the search for the nuclear-dibaryon phase- 
transition at densities ranged from 1.5po to 4po- This requirement constrains the variation 
of dibaryon mass to 1940MeV^ < Md < 2mQMeV. 

A comparison with the hadronic field theoretical model of reference shows that 
in the Hartree approximation a nuclear-dibaryon heterophase appears rather than a pure 
condensate. The Bose-condensate for d' dibaryon appears at 3 times the saturation density 
of nuclear matter. Independent calculations using the quantum field theory of hadrons 



and dibaryons have been developed in In this framework the transition density for 
nuclear and neutron matter is Pd/po = 2.87 and Pd/po = 2.16 respectively ||22|, similar 
results are found in reference po/Po = 2.69 and p/j/po = 2.57 for nuclear and neutron 
matter, respectively. Therefore our prediction of Pd/Po = 2.55 is consistent with these 
values. 

In this work we have studied the properties of a system composed of symmetric 
infinite nuclear matter with dibaryons. This has been performed in a theoretical frame- 
work which takes into account the finite size as well as the quark structure of nucleons 
and dibaryons. The variation of particle properties with baryon density have also been 
considered. The nucleon mass shows a monotonic decreasing behavior for all the densities 
and dibaryon abundances Q studied here. At Q = nucleon swelling is observed at den- 
sities below the normal nuclear matter density. The bag radius increment is less than 2%, 
in accordance with theoretical estimates based on y-scaling interpretation of quasi-elastic 
electron scattering A nuclear-dibaryon matter phase-transition is found at zero tem- 
perature. The dibaryon condensate predicted in this work could be an intermediate state 
before the transformation to quark matter. If a dibaryon condensation takes place for 
baryon densities in the range 1.5 < psl Po < 4, then values 1940MeV^ < M^) < 2000MeV" 
are predicted for non-strange dibaryon masses. As a consequence of phase transition the 
compressibility of the system is considerably reduced at high densities, as it should be 
expected from astrophysical scenarios. 

Extensions of the present work in order to describe strange dibaryons can be in- 
mediately performed by inclusion of the color electric and color magnetic interactions in 
the MIT bag model, also multiquark bags describing some exotic states of multibaryons 
could be treated in this framework. 
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Mi) [GoV] pilfiii pnlfiii hh/ki 



1.94 


0.228 


0.335 


1.30 


1.95 


0.283 


0.380 


1.09 


1.96 


0.335 


0.429 


1.01 


1.97 


0.383 


0.477 


0.98 


1.98 


0.427 


0.522 


0.96 


1.99 


0.470 


0.566 


0.94 


2.00 


0.505 


0.603 


0.93 



Table 1: The baryon density p and the quotient of thermodinamical compressibihty 
K — dP/dpB evaluated at the transition point as functions of the in- vacuum dibaryon 
mass Md- The label I (II) corresponds to pure nuclear (dibaryon) matter state. 



Figure 1: B (dashed hne) and the nucleon bag radius, for the values Q — 0, 0.25, 
0.33, 0.47 and 0.5 (solid lines) as functions of the baryon density, pb- The arrow indicates 
the increasing Q values. is the nucleon radius at zero baryon density. 



Figure 2: The dibaryon chemical potential (dashed line) evaluated aX Q = 0.5 and 
two times the nucleon chemical potential (solid line) at Q = plotted as functions of 
the total pressure of the system. The phase-transition for in-vacuum dibaryon mass 
Md = 1970 MeV occurs at a baryon density of p/ = 2.55 fm~^ in the nuclear matter 
phase, corresponding to a baryon density of pn = 3.17 fm~^ in the pure dibaryon phase. 



Figure 3: Pressure as a function of baryon density for the dibaryon abundances Q = 0, 
0.25, 0.33, 0.47 and 0.5 (dashed hues). The arrow shows the increasing Q values. Sohd 
line represents the physical states of the system, corresponding to nuclear matter (Q — 0) 
below the relative density ps/Po = 2.55 and to pure dibaryon matter {Q = 0.5) above 
the value Pb/po — 3.17. The horizontal segment between these densities represents a 
coexistence region between the two phases. 
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